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Accurate ocean thermal structure information is an impor- 
tant requirement necessary to support U.S. Navy anti-submarine 
warfare (ASW) operations. Estimates of oceanic conditions 
that are to be used to generate sound velocity profiles and 
acoustic range predictions can be obtained from climatology or 
objective analyses. These sources of information are usually 
represented digitally on coarse grids which limit their 
spatial resolution. Neither climatologies nor objective 
analyses can depict the true oceanic conditions. 

The best possible source of lecal information for a ship 
at sea is an observation. The oceanic conditions at a remote 
point, however, are often essential for mission planning. 

The question that arises at sea is how can a local observa- 
tion be used to estimate conditions at some distant point. 
In this thesis, a simple technique for accomplishing this 


will be developed and tested. 


A. EXTRAPOLATING A LOCAL OBSERVATION TO A REMOTE POINT 
Some attempts have been made operationally to “enhance" 
a climatological temperature profile at a remote point by 
adding the departure from climatology determined from a 
local observation to the climatology at a remote point. 
Such a technique has been reported to provide useful esti- 


mates for computing acoustic ranges in Navy fleet exercises 


(LT D. Pedneau, personal communication). This type of 
extrapolation is also used eenesispanded Oceans inecrmrat 
Structure analysis (EOTS) “Holl, et val. 1972 Seetpe ed a - 
Fleet Numerical Oceanography Center (FNOC). A preliminary 
step in EOTS is to estimate values at grid points from near- 
by observations using this method. 

Various methods of extrapolation also constitute the 
fundamental approach used in other objective analyses. 
Simple methods such as used by Druyan (1972) involve assigning 
values to grid points by weighting the observation inversely 
proportional to the square of the distance from the grid 
point. Other more sophisticated methods require the use of 
a "trial value" at any arbitrary location (Kruger, 1969). 
The concept ofa "trial value" is used in the same context 
that a climatology is used in this thesis. In these applica- 
tions, the deviation of the trial value from the observation 
(climatological anomaly as used here) is multiplied by a 
weighting factor before adding to the trial value (clima- 
tology) at the extrapolated position: oprecaine analyses of 
this type differ according to the scheme used in determining 
the weighting factor. One wide1ty used scheme, optimal inter- 
polation (OI), obtains a weight function through autocorrela- 
tion techniques (Alaka and Elvander, 1972). OI is used 
widely for meteorological analysis (Bergman, 1979; McPherson, 
et al., 1979) and is being adapted for ocean thermal analysis 


applications at FNOC (Innis and Williams, 1983). 


A Navy ship at sea may not have a large number of 
observations available to generate a regular lattice of 
representative values for its operational area. Thus, util- 
izing a single observation to estimate oceanic conditions at 
distant locations is a desirable approach. The simplistic 
method of enhancing climatology described earlier is really 
a special case of OI when only one observation is used. [In 
this thesis, a method of extrapolating an observation in 
real-time to a remote position will be studied. This method 


will be referred to as "simple enhancement." 


Bem OCLANTC TIME AND SPACE SCALES 

The accuracy of ocean thermal structure estimates 
obtained using simple enhancement depends a great deal on 
homogeneity of the oceanic region under consideration: 
Knowledgeable use of the observations thus requires under- 
standing of temporal and spatial scales of oceanic 
wat taal 1 tye 

Significant ocean thermal structure anomalies can range 
in size from 100 km to the size of the ocean basin (TOPEX 
pelemece Working Group, I98l). For instance, in the North 
Pacific, large oceanic and atmospheric anomalies have hori- 
Zoe Semiemsmens the Order of 1/3 to 1/2 the size of the 
basin (Namias, 1972). More recent studies by White and 
Bernstein (1979) used autocorrelation analysis to determine 
the space and time scales of Vomrabmieltyein the North Pacific. 
Zonal length scales were found to be 1500 km with the 
meridional scales approximately half that value. The zonal 
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length scales decrease to about 300 km at 300 m depth in the 
western portion of the mid-lattende Novenee ene eee eres 
indicates greater variability in the western region. 
Meridional wave number spectra of temperature from hydro- 
graphic sections in this region (White and Meyers, 1982) show 
similar results of 400 to 750 km length scaleswat depths of 
100 to 300 m. 

These statistical results give only a crude indication 
of the horizontal space scales for which simple enhancement 
may be useful in ASW applications, i.e., several hundred 
kilometers or more. The existence of mesoscale eddies in the 
North Pacific (Bernstein, 1974; and Bernstein et al, 1982) 
will complicate the method. Large errors would be expected 
when extrapolating’ an observation across the boundary of such 
an anomaly to obtain an enhanced temperature profile. Further 
sources of error are realized as a result of relatively short 
time scale changes in the vertical structure. 

Urick (1975) noted the velocity of sound near the surface 
1s sensitive to local changes in the temperature profile due 
to heating, cooling and wind mixing. Watt and Morrice (1980) 
discussed the tactical significance of sea state and near- 
surface temperature and salinity structure and noted its 
dependency on atmospheric conditions. They also demonstrated 
the significant difference in sound ray paths between morning 
and afternoon under condittiens of TPidhtwwainds sineesereug 


Surface warming. Intense surface thermal stratification due to 
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diurnal solar heating is also documented by Shonting (1964), 
who found variations in sound velocity of 2 m per second in 
the upper 5 m of the water column. 

Investigations on upper ocean response to wind mixing 
(Elsberry and Raney, 1978) show variations in mixed layer 
depth (MLD) and sea surface temperature (SST) with time 
scales of several days. Similar studies by Elsberry and Camp 
(1978) show large changes in SST are normally accompanied by 
changes in MLD and occur in association with periods of 
SerongmatLmnospherme LOorcing. Events Of these types can 
produce significant departures from climatology (Camp and 
Pas berry, 71973) - 

The feasibility of simple enhancement can now be more 
clearly examined. Based on the preceding discussion, the 
method may seem incapable of demonstrating any skill in 
accurately depicting a temperature profile. The real skill 
may become more apparent in depicting a "representative" 
temperature profile for a given region. This may be more 
appropriate for ASW applications. It is of significance that 
active sonar systems with typical ranges of 30-40 km make an 
implicit assumption that the ocean is homogeneous over these 
ranges. Thus, range-independent acoustic models usually 
require a sound velocity profile at a single point. Due to 
internal waves and other small scale fluctuations, an 
observation of a temperature profile at a single point may not. 


be representative. The question becomes, "Is an accurate 


ame 


depiction of MLD critical OF would an average Er eee 
region be more suitable?". Here the point can be made that 
the usefulness of simple enhancement will depend on the type 
of application for which 1€ will be used>) (22 ome Bs mace 
observation is available, then an average MLD cannot be 
obtained. As additional observations are taken, however, 
increasingly more sophisticated methods can be utilized to 
obtain estimates of the ocean thermal structure. The experi- 
ments in this thesis will focus on the use of a single obser- 
vation. The last experiment will address the feasibility of 
obtaining a "representative" temperature profile by using a 


filtering technique described in Appendix A. 


a? 


Ii. FORMULATION OF SIMPLE ENHANCEMENT 


A. USE OF A SINGLE OBSERVATION 
Peony lo Obsernvaceem,s imi!) , at position lecan be used 
to enhance the climatology, C(x) at some other location, x, 


by the following formula: 
Ce(x) - C(x) + [Tm(1) - C(1)] (1) 


where Ce is the enhanced climatology and Tm consists of 


Speenvaumlonaleerrors, ©, plus the true value, T(1l). That is, 
Tm(1) = T(1) + e€. An alternate relation is: 
Geos) = tm) 4 Ge, sc = Gx). — C(1). 2) 


Here it is assumed that a suitable climatology is available 
and a value at any location, x, can be interpolated from the 
Shinacolody-— oluse, a z2eal-time ocean thermal structure 
analysis or forecast (Clancy and Pollak, 1983) could just as 
easily be substituted for climatology as the "trial value." 
Equation (1) provides a convenient means of enhancing 
climatology using an observed anomaly. The equivalent form 
shown in equation (2) indicates the dependence of simple 
enhancement on the trial value used. For example, compare 
the use of a synoptic forecast with that of climatology as 
trial values. Two cases can be examined. Case (1): The 
Forecast and climatology fields may be dissimilar in magnitude, 


but show the same general trends of large horizontal gradients. 
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That is, one field may show a nearly consistent bias over the 


other. In this case, the 6@eterm woulesbewsamn ciao see) 
the climatology and forecast.) Enus tne) values wre rsrcen 
would show little difference. Case (2): The forecast field 


is represented on a grid finer than the climatology, and 
provides greater detail than the coarser grid. The 6C term 
would then be different for both fields, and the enhanced 


climatology and enhanced forecast would be different. 


B. ESTIMATING THE ERRORS USING GCP PiMUM STE Reon Leh 
The extrapolation error, I which results from simple 


enhancement 1s: 


I, = T(x) - Ce(x), (3) 


where T(x) is the actual temperature. This relation will 
prove useful in experiments testing the feasibility of simple 
enhancement; however, a means of theoretically estimating 

the expected errors prior to performing experiments would be 
desirable. Thus’ the theory of optimum interpolation which 
minimizes the mean square interpolation error is examined. 


The basic equations (Alaka and Elvander, 1972) are: 


[Tim (je (a) EG see (4) 


@ 
Team n Cea) ie 
j=l 


where Io is the interpolation error andmie the numecame: 
observations. The weights, P(j), give the relative importance 


of each measurement, and are determined by: 
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nN 
MPVEVE TT Ge meme) = yp (x,i) i = 1,2..-n (5) 
j=l 


and the minimum mean square error, E, is: 


5 Nn 
E =o [1- £ u(x,i) P(i)] (6) 


i=l 
where 6? is the signal variance, 1/A is the signal to noise 
ratio, and u(x,i) is the autocorrelation between T(x) and 
eos tinalivy wid, Jers the autocorrelation between T(1) 


meni wiichmamvOolyes pDroaucts Of all pairs of observations. 


poe N=) equations {4)—(6) reduce to: 
moo) —seGom- (im(l) = C(1)] Pp + Io (7) 
P= ny ee). where u(l,l) = l (8) 
and 
E = o°(1-p(x,1) PI (9) 


with subscripts for P and »} dropped. Setting P = 1 in (7) and 
comparing to (1) shows that simple enhancement is a special 
case of optimum interpolation. From (8), P = 1 implies 
u(x,1l)=1 and \S = 0. The assumption is that the correlation 
between point l and point x is perfect and the error of the 
observation is zero. For simple enhancement, OI provides no 
means of estimating the errors since (9) reduces to the 

trivial solution of zero. If a value other than one is calcu- 
Pree eter yee sme tnod can noe Longer be called Smo Le sonhence— 
ment. This variation of simple enhancement with P=1l and n=l] 


Wile be called “OL enhencement." 


ile) 
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A. GENERAL APPROACH 

The concept of simple enhancement was tested by using two 
different sources of synoptic fields for trial values. 
Expendable bathythermograph (XBT) data were used to provide 


"locally observed" temperature profiles to be extrapolated to 


the enhanced position. The data set also provided a control 
for comparisons against enhanced values. Equation (1) was 
applied and the results evaluated with equation (3). To test 


dependence on the trial value, both a climatology and a synop- 
tic ocean thermal forecast were used in the experiments. 
l. The Data 

The data used in this study were provided by the Naval 
Oceanographic Office. They consisted of a subset of XBT drops 
made by the USNS Silas Bent while surveying from Kodiak, 
Alaska to Hawaii. The sampling started on 26 September 1982 
at 54°N, 149° 30° W, continued on) an tbackm couche ee 
41° 45'N, and ended 1 October 1982. Observations were made 
approximately every 15 km using 750 m Sippican XBTs with 
temperature and depth accuracies of 0.2°C and 1 percent, 
respectively. Data were recorded on magnetic tape and analog 
recorders. 

Initial reduction of the data was done at the Naval 
Oceanographic Office giving temperatures at 1 m depth inter- 


vals. For this study, additional screening of the data was 
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made by plotting temperature profiles of the 85 XBT drops 
made. Only one XBT observation was removed from the set based 
on what appeared to be "wire stretch." This is a condition 

in which the unravelling wire of the instrument is stretched, 
causing increased resistance in the wire and spuriously high 
temperatures. 

A grid was constructed that divided the north-south 
track into 100 equal spaces and the data were linearly inter- 
polated to it. The interpolated reports are numbered 1-100 
from north to south for reference purposes. 

The resulting vertical temperature cross-section from 
the surface to 300 mis shown in Figure 1. Large horizontal 
and vertical variability is apparent in the upper portion of 
Pioesedectmalettermogline near 50 m. This variabi¥aity is 
presumably due to shallow mesoscale eddies, fronts and inter- 
nal waves. The domain includes most of the subarctic transi- 
tion zone which has been studied extensively by Roden (1970, 
1971, 1977). The colder water below 100 mis typical of the 
intermediate water mass of the Pacific (Reid, 1965). A 
common characteristic of the data is the appearance of temper- 
ature inversions between 100 m and 150 m (Figure 1). These 
inversions agree with observations by Roden (1977) in nearly 
the same region. 

Mgt eoecaeTOanwOr tne Subarctiec front 1S farely persist- 
ent on a weekly time scale (Roden, JES) diet ivoreal horrzoncal 


temperature gradients in the frontal region are 1-2°C/100 km. 
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Two such fronts with gradients of 1.9°C/100 km occur between 
about 45° 20' N to 45° 45' N and between 43° 10' to 43° 45' N 
(Figure l). 

Above and below the thermocline there is almost no 
vertical coherence in the thermal field. Except for the 
Variations in the 40 m and 60 m region, the horizontal length 
scale is not easily determined by inspection. However, the 
scale below 60 m can be estimated to be at least the length 
of the survey track (about 1200 km). This is consistent with 
length scales in the eastern North Pacific described earlier. 
In all respects, the data set is typical for this season and 
region of the ocean and should be sufficient for a feasi- 
bility test of simple enhancement. 

2. The Thermal Fields 

The thermal fields used to provide trial values in 
this study were a monthly climatology and a daily forecast 
of the northern hemisphere upper ocean thermal structure. 
Both fields are represented on the northern hemisphere 63 x 
63 polar stereographic projection used at FNOC. 

Climatological temperature profiles were extracted 
from the monthly climatology used at FNOC. The climatology 
was constructed from the surface to 400 m using data from 
the National Oceanographic Data Center (NODC), and data 
received through direct exchange from foreign Navies and 
foreign and domestic institutions (Bauer, 1982). Its struc- 


ture 1s represented vertically on a fixed level grid. The 
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top of the permanent seasonal thermocline and its shape are 
described separately by a set of variable levels (Table l). 

Forecast temperature profiles were extracted from 
fields produced by the Thermodymamic OcéanyPr edie en syseen 
(TOPS) run daily at FNOC. TOPS cycles with@a veal eim- seca 
thermal structure objective analysis in an analysis-forecast- 
analysis fashion (Clancy and Pollak? 963) —and@waS Soece@e--c 
and evaluation phase during the period of this study. The 
verticalMgrid is shown 1n»ytaolewms 

The TOPS forecast presently uses the Level-2 turbu- 
lence closure of Mellor and Yamada (1974) to parameterize 
the vertical diffusion of momentum, salt and heat within the 
mixed layer. The parameter of this model that controls the 
MLD evolution is the Richardson number. The base of the 
mixed layer is the depth at which the local Richardson number 
exceeds a critical value. The Level-2 closure differs from 
those used in a bulk model, such as Garwood (1977), where 
MED 1S “explicr ely (pe ecdteced- 

Vertical temperature cross-sections from the surface 
to 300 m (Figures 2-3) are shown for the climatology and 
forecast fields. 

Monthly climatologies for September and October were 
extracted and linear interpolation was used to obtain values 
for 29 September, the center date of the data set. Clima- 
tological profiles were then interpolated to the 100 equally 


spaced data points. This interpolation was nonlinear using 
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Table l. 


a. 


Vertical grid used to represent the TOPS forecast 
temperature profiles. 


Vertical grid used to represent the FNOC climatology. 
In addition to fixed levels, a set of floating levels 
is used to represent the location and shape of the 
top of the seasonal thermocline in the climatology. 


a. TOPS b. Climatology 
Depth (M) Depth (M) Floating 
+) Oe OF Primary layer depth (PLD) 
25 252 Temperature at PLD 
Ue 7 OE, Temperature 25M above PLD 
2 ion eS 12.5M below PLD 
T25 ILO tir oie 25.M below PLD 
Zou 25.6 Ge 50.M below PLD 
32.25 ESO 
40 200 
50 250. 
6225 300. 
ioe 
eZ 5 
ales) 
200= 
30G2 


*Not a computational level in TOPS numerical integration. 
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a standard subroutine (INTRPS) residing in an FNOC subroutine 
library. Finally, linear interpolation was@useaqsvereuca 

to obtain values at 1 m intervals to match the data. The 
climatological profiles were then subtracted from the observed 
profiles to obtain anomaly profiles. A similar procedure was 
used to compute anomaly profiles for the TOPS forecast. 

Since the effect of diurnal and other short time scale changes 
will not be studied, a single TOPS 24-h forecast will be used 
for the entire 7 day period. During the cruise, XBT observa- 
tions were sent to FNOC and assimilated into the daily 
analysis/forecast system. Comparisons with data are nonthe- 
less independent since the TOPS forecast is valid on 26 
September, the start of the data set. 

Figures 4 and 5 show contours of the climatological 
and forecast anomalies. Note that the observed temperature 
field differs from climatology by several degrees C. Such 
anomalies are not uncommon. In the forecast, the differences 
can be explained by imperfect initial conditions, anaccurate 
atmospheric forcing, or inadequate parameterization in the 
mixed layer forecast. Regardless of the reasons for these 
deviations, their existence suggests that a technique to 
reduce the size of the anomalies, such as simple enhancement, 


would prove useful. 
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Same as Figure 4 but for TOPS forecast anomaly. 
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B. ENHANCEMENT USING A SINGLE OBSERVATION 

The 100 data points were subdivided into three over- 
lapping sections for detailed study; points 1-40; 31-70; and 
61-100. To test simple enhancement, the anomaly profile 
(observed minus climatology) at the first point of each 
section was added to all the climatological profiles in that 
section. The first profile for each section exactly repro- 
duces the observation and adjacent altered profiles are called 
"enhanced climatology." A similar procedure was used to pro- 
duce "enhanced TOPS," except that the forecast anomaly 
(observed minus forecast) was added to all forecast profiles. 

Error fields were computed (observed minus enhanced) for 
enhanced climatology and enhanced TOPS and contoured from 
the surface to 100 m (Figures 6-11). Large differences 
occurred between 30 m and 60 m. This indicates that the 
vertical temperature gradient at the base of the mixed layer 
cannot be accurately extrapolated forward in space using this 
method. The fact that enhanced temperature profiles for 
climatology and TOPS are nearly the same iS more readily 
apparent in examination of individual profiles, which are 
Shown for every fifth point in Figures 12-14. The similari- 
ties of enhanced climatology and enhanced TOPS is consistent 
with equation (2) and the discussion of Case 1 in Section II. 

Notice also that erroneous temperature inversions were 
Ptensedncedmuer section l at™points 6, 119P 16, "21, and 26 


(Figure 12). The reason for this can be explained by an 


2) 


‘doa, 4B pezeotpUtT 
(T UOTIOSS) OP-T SRzUTOd TOF Anq ‘Gg pue FY Sainbty ut sattTewoue se Aem oues 
94} peanojzUuO. pue pejnduods AboToOReUT[TO padueyuse Jo pTety AoA °*9 JINbHTYy 


SalI] Lo 





(WIHL dd 


28 


POMUOFeso- PUEM/—[¢€ SjzUTOd 


Lv 


TOF 3nq g9 sANbTty se sues 


SIC ee) 9) 


QV 


GP 


“f oanbty 





208 22d 


882 


S 
& 


86 


AS 


BZ 


a9 


8S 


av 


ihe 


82 


ol 


(WIHLdS0 


ay 





30 


42 


43 


44 


45 


Be | Ube 
same as Figure 6 but for points 61-100 and section 3. 


Figure 8. 


"Jsed0103 SdOL pooueyuea AOFJ Anq 9 SANbTy se swes °*6 OANbTyY 


JUNL TL | 


WS 1S 25 ES 


\ 





Jao 





AG 


Ag 





0/ 
a9 
as 
ar 
pe 


( ee 
os az 


Al 


N 


z2a 


' = ros) 


HP “ts Hh? Ht 


a 


(WIHId 


oa 


S 
| Ag Aa A LA 
C™ a4 WTS 





082 2 
Bi] 
69) 
O82 © 
O 
08a ) 
My 
(2) 
uy 
4y"1%) 
3 Ww 
O 
= > 
re) 
oe) 
O 
o 
40) 
& 
S 
0) 
3 
Lig 
O SH 
Bo) } 
S eae = 
ne ae 
— 
2s m 
225 =; q) 
My 
= 
'e)) 
To 
fy, 
6 
= © 
2 
0) 
= 
7 10) 
_ ep) 
282 
ier, 7 
<E 
at 
co) 
\y 
a 
'e)) 
or 
fx, 
02 Cge@ Jag 22 ao 
S ~ 3 3 3 3 3 2 3 2 3 


hua 


= a - ~ WES 4 26, ec — = pi = o 
} f 
ue —_ x ges ee on ere Ut~. Coat ae 4 ee La eu 
ha ~s ~ we eee —— =, ~, ~ — -- 
ad a = Oo oo" ~ ~e NO ee ee — fj 
reese ee) *. oo 
— . 4 ~ nen a 
f = ate sy . —_* eee ne 
e a *. —_ eo” 
- een = > 
- wom . “ 
- ~ 
f *¢ 
| - s- eS 
' “a1? 
ate at ee 
ee 
{ef ‘ ——— 
-1t Sey a ee ee, 
_. . Niet - ae ere 
. _ 
als ~ oo 
_ — ~ 
L oe 5 “cm. . ‘ mS) 
- - AL *. o +s 
f \ US ~ —_—-" 
‘ . 
‘ —_ ‘ 
‘ ——e 1 
i - —s \ 
/ ‘ —— ~~, \ 
j ' ' oie \ 
/ 
/ 
' 
4 aa 
\ *, 
: ; eae at 2 . aoe 
f -———s, o~ 


. 
a pe ~ ts 
. ~ «. 
“= ———* . —_—— ae * - 
Me -~ ~ — ~ —_ e 
re a ™ — 2 om 
oe <s., “ea tlel, =. . = aaq ' 
——— = ~~ te [a , ———— ~~ & bar 
—_— ~. 5 e —= 
- “ = ae 








cy) a ‘ | — — 
rs a 
s / \ a — 
r -—= \ ———— 
4 —— 
j / _* oe 
4:2) Soto — Ee P X 1 \ - — ' 
rai — es as 2 eae ——— Cai ake a3 + 
a pee Fearne cee ee — “wa k ‘ ! 
ea a * ees eso_—: } ‘ } — 
== _ i 00s cere 3 ‘ - —_—_—_—_—_ 
aie | ae Rac F 
jc tn AV, wom et a 
i = [ yh , —, \ 4 t 3 
Pettitt? ome pa esN 
Pe c t i Pe \ o ‘ tee, 
= ) SU \ }\! 1 ae —— { ~.. 
a | \ rij 3 = aaa So as 
' \ ~ Hewat - an if rs : 
, f ” e a atl 
oe st ae 
pay 
Nee eC Che. id oval . Wy <- mike 
at = 
ee ee +. ~~ eee & ——s., a Ne ean, a” 
er ——— a —aly oe rae 
we ‘ a .’ y bat LS f - ‘~, 
- Sree ty . : PN 
- a Se eerie Ns 
. of 4 oy | wo? _” - ——_ So to 
~ a Pe re . _ o* * . a 
ee ee — “ oe —— _ . 
=-3C a -= —a-a en a \ 2 
wei ax = A= : Fe, 34 
: BE ———— a5 spenae! G eg 
= 
1 < ‘ 
“vals Fl ' jz. on a eS <o 
os 1 Fis — — 
1 Ca | ° —_—_— 
=—* —_—_ a 
1. ee = = 
ao / joe Fi st ol 
“Sf _ 1 Se ota 7 
~o-L7 oon 
a! \ 
o \ 
- 
Choe DAG, 
; ae (| 
“+ ae me 
, Bs 
—_—_——e a ib 
~ — ! 
Oe ae hm “ ~ . 

-- oe “. = ( 
= 7 | ls *~ Fe 
zt tif wm, ’ | 
: = ba 0) 

| | : He ' * 3 Se Ae p Tad on 

“. Pid > = 

- ' — a qs 
SE me Saye a eam 
SS a . * £ — 
- 
ice eae A = 
a epee ae it yy ‘ 7 ' ~ a. : 
- a — 
“e a7) a ‘ 





ve 
—_ ce | _* at - ase 
z 265 _—--- A 
eae f eee 
7 el — 
| ‘ | ss =“ —— oe 
} s f ! - ot me a 
| WE ope { ; vie cass ere 
' f ~~ 
é Wow! J a ace 
-—? 
| { ‘ t \ oa 
4 | { i / 
| re 
e —. 


dW ga . Add 


) 

) 

" 
10uU 


a>) ia) ic? i= 
Caan] oar, 


60 
8 
2 


€ 
Us 


(WIHLs su 


33 


LHI TP LIbt 


Same as Figure 8 but for enhanced TOPS forecast. 


Figure ll. 
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example showing the steps of simple enhancement. First note 
the depth of the inversion at) point 11 (Figure 12) a. 
indicating this depth is shown at points 1 and 11 through A, 
B, D, and E. The difference between the observation and clima- 
tology at this depth is about 2°C. Adding Ehis dieeerences--= 
D gives the enhanced temperature E at point 11. Compare this 
to the enhanced temperature G computed the same way only 

using C, B'and F at a different depth. The apparent inversion 
due to the higher temperature at E results from the fact that 
small differences 1n MLD may be associated with large differ- 
ences in temperature in the MLD region. When these large 
temperature differences are added to an independent clima- 
tological temperature profile at a different location, the 
possibility of an inversion results. 

For each section, the mean and RMS errors (Tables 2-4) 
were computed for depths of 0 m, so0em, )  Olmemerencls)o Omen 
Confidence limits were calculated using a student T distribu- 
tion with n=40. In all cases, theverror 16) certineaeas ene 
Observed value minus the trial value. To composite the 
results, averages for all three sections over depth were com- 
puted (Table 5). The mean errorsfor enhanced climatology and 
enhanced TOPS were nearly an order of magnitude smaller than 
those for unenhanced climatology. The average RMS errors for 
the enhanced values were about 1/2 that for unenhanced clima- 
tology. Also, the errors in the enhanced climatology are 


about the same magnitude as errors in the enhanced forecast. 
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Thus, two disSimilar trial values produce similar enhanced 
values with nearly the same improvements. This is a desirable 
resule Since the choice Of the trial value used was not 
emaseica lL. 

Since simple enhancement was shown to be a special case 
of optimum interpolation, the former might be improved by 
including the known spatial statistics of ocean thermal vari- 
ability via this technique. To test this hypothesis, equations 
(7) and (8) were used to compute an OI enhanced climatology. 
The signal to noise ratio, 1/A, and the autocorrelation, 
u(x,1), could be computed from the 84 observaticns. A more 
desirable approach would be to obtain the constants from an 
independent and much larger data set. From White and Bernstein 
1 ge9) , ae ~oegepeey4 = .432. fhe autocorrelation was deter- 


moneda from 


2 
oo ee q 7c fAX) (10) 


(See Clancy, 1983), where AX is the latitudinal separation 
between the observed and enhanced value. The constant c was 
determined by fitting (10) to the autocorrelation curve in 
White and Bernstein (1979). Although these statistics are 
depth dependent, values for the sea surface were used and 
Sooltee touabimceptns tor simplicity. Table 6 lists the 
resulting P'S as a function of AX. Equation (9) was used to 
compute the expected minimum error for each enhanced profile. 
The average minimum error for all profilles (Tebie 5) 1s OFO2°C, 


compared to the actual mean error of 0.80°C. 
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Table 2. Mean and RMS temperature error (observation minus 
test value) for test values of climatology, enhanced clima- 
tology, enhanced TOPS and OI enhanced climatology at depths 
of 0, 50, 100 and 150%™ for section | (pores {7 


Mean Error °C With RMS Error °C With 
95% Confidence 95% Confidence 
Intervals Intervals 


Deets Orin 

CLIMATOLOGY ee! og osy TOC oe pee) 
ENHANCED CLIMATOLOGY -.42(-.52,-. 32) {52% 2427262) 
ENHANCED TOPS meee OS) ONS) .90(.40,.60) 
OI ENHANCED CLIM. OL (= ese a Gane. | Gaus 35) 





DEPTH: 50 m 

CLIMATOLOGY ee OO ely OT reo) TIO (TTS , 2e24) 
ENHANCED CLIMATOLOGY =a (202 Pay a, ey ARP 25651610.) 
ENHANCED TOPS -1.58(—2.11,-1.04) 2.30(1. 76 722463) 
OI ENHANCED CLIM. =) -18(-1. 737 =565) 2.09(1.54,2.64) 


DEPTH: 100 m 
CLIMATOLOGY a AE eS) 


ENHANCED CLIMATOLOGY 03(=,037.03) 
ENHANCED ‘TOPS oO 05 7 ks) 
OI ENHANCED CLIM. -33(=e0,-- 26) 


DEPTH: 150 m 

CLIMATOLOGY .54 (-.59,-.49) 
ENHANCED CLIMATOLOGY 1S 13).20) 
ENHANCED TOPS 33 (2307220) 
OI ENHANCED CLIM. » | -.14(-.21,-.07) 
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Table 3. Same as Table 2 except for points,.31=-70 and 
SEGeVON = . 


Vecnmeeaaneor -C With RMS Error °C With 
955 Conti aence 95% Conti dence 
Intervals Intervals 


DEPTH: Om 

CLIMATOLOGY -44(.40,. 48) -46(.42,.50) 
ENHANCED CLIMATOLOGY AO ees) e202 0) 
ENHANCED TOPS - LO { eee 2 4) 225(.203ae1) 
OI ENHANCED CLIM. w20(.22 pee) 126 Wee OL) 





Deen: 50 m 

CLIMATOLOGY .11(-1.66,-.57) -04(1.4972859) 
ENHANCED CLIMATOLOGY (OA reea06) ietdeZ4, 2.33) 
ENHANCED TOPS .49(-.04,1.02) TA (2052. 298 
OI ENHANCED CLIM. =a cOn. 37) . 84 (1.26 ,2. 48) 





DEPTH: 100 m 
CLIMATGrROGY -1.03(-1.15,-. 90) lL Gee pelle ae 
ENHANCED CLIMATOLOGY = .0/7(=.20;505) HOO. 27 852) 
ENHANCED TOPS = .02(=3137 790) -50(.249eey) 
OI ENHANCED CLIM. S04) OO ee ogeoS) 





DEPiInH: 150’ m 
CLIMATOLOGY ~10(-1.25,-.94) i 2004) 1. 36) 


ENHANCED CLIMATOLOGY ~26(—=42,-. 10) 50 40../ 2) 
ENHANCED TOPS eS se OL ee -48(.33,.63) 
OI ENHANCED CLIM. 64 (=e 5) .86(.68,1.05) 
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Table 4. Same as Table 2 except for points 61-100 and 


SeCELOnmo. 


Mean Err@r °C With RMS@bisewece with 
95% Confidence 95% Confidence 
Intervals Intervals 


DEPTH: Om 
CLIMATOLOGY 

ENHANCED CLIMATOLOGY 
ENHANCED TOPS 

OI ENHANCED CLIM. 


DEPih <= 50 m 
CLIMATOLOGY 
ENHANCED CLIMATOLOGY 
ENHANCED TOPS 

OI ENHANCED CLIM. 


DEPTH: 100-m 
CLIMATOLOGY 
ENHANCED CLIMATOLOGY 
ENHANCED TOPS 

OI ENHANCED CLIM. 


BEET el 50 ein 
CLIMATOLOGY 
ENHANCED CLIMATOLOGY 
ENHANCED TOPS 

OI ENHANCED CLIM. 


56 (.48, .64) 
-34(.26,.41) 


eo US 23} 3) 
43 (..35 7252) 





-39(=2,95,=225) 
1 OSS) 
-14(-.06, .34) 
Rey Gaul Ik pe 7S), 


‘26 (= 157, eee) 
-3O(e20y ee) 

. 30( 22041) 

39 (= aera 4) 


.41 (=1 252730) 


Sear ce Clan 
~19(=.29 aoe) 


-69 (=. 79 jae) 
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~Ol. (253 72iGa) 
42 (.34,.49) 
.36(.26, .46) 
014 342,,.59) 


2.43(2.23, 2568) 
68 (944, 68) 
-63(.44 res) 
sien oT re ooy 


- 30 (1. 20a) 
polaaey 02) 
-44(.34,.55) 
“fOr oO pea) 


749:(1.. 3474556) 
B30 (- 2am) 
-38( . 2iqmee Sy 
-/Ol 66 poo) 





Table 5. Average mean and RMS error for all three sections 
and depths. 


Meal) Error RMS Error 
(2c) (Pe) VE 
CLIMATOLOGY Tae) 
ENHANCED CLIMATOLOGY 267 
ENHANCED TOPS 66 
of ENT ANGE De eiM. 7 OO Oe 
FILTERED ENHANCED CLIM. ee (0 
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Table 6. Weights, P, as a function of latitudinal 
separation, AX, used to Compure an Ol enhanced climatology. 


AX (Deg) P AX (Deg) P 
0.0 0.6981 Desi 0.6040 
0.1 0.6979 Da 0.5955 
OeZ 0.6972 Des 0.5868 
0.3 0.6961 2.4 0.5778 
0.4 0.6945 2S 0.5686 
0.5 0.6924 2.6 0.5591 
0.6 0.6899 Dow 0.5495 
One 0.6870 Ie 0.5396 
0.28 0.6836 2.9 0.5296 
0.9 0.6798 30 0.5194 
1.0 0.6756 2a ik 0.5091 
ivan 0.6709 3.2 0.4987 
122 0.6659 Be 0.4882 
1.3 0.6604 3.4 0.4776 
i, @! 0.6546 3.5 0.4669 
1.5 0.6484 3.6 0.4561 
1.6 0.6418 cae 0.4453 
1 0.6349 3.8 0.4345 
8 0.6276 3.9 0.4236 
1.9 0.6201 4.0 0.4128 
20 0.6122 
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THe mean and RMS ereors (Tables 2-5) indicate that no 
improvement is realized by including spatial statistics in 
simple enhancement. This is probably due to the use of a 
Single observation for enhancement. Although the expected 
error is reasonably close to the actual error, a greater 
improvement might result if multiple observations were used 
fimcombination with equations (4) and (5). In addition, it 
should be appreciated that the OI technigue guarantees 
Minimal errors only in an ensemble-mean sense. Thus, one 
particular realization, such as done here, might not neces- 
Sarily show the advantage afforded by this technique. 

Gee OS AINING A’ REPRESENTATIVE ESTIMATE OF THE OCEAN 

THERMAL STRUCTURE 

The fluctuations at the base of the mixed layer noted 
earlier can be removed by a low-pass filter while retaining 
large scale trends of a few hundred kilometers or more. A 
Simple procedure to "correct" each temperature profile 
using a 5-point running mean horizontally (see Appendix A) 
was applied to the data for the three sections separately. 
Anomalies from the filtered data were then used to compute 
enhanced profiles for climatology as before. Since the 
mixed layer depth was averaged over 5 adjacent temperature 
profiles, small scale fluctuations will not be present in the 
filtered data. The resulting enhanced temperature profiles 
would then be "representative" over the 5 points, or about 


48 km. The contoured error fields of enhanced climatology 


47 


resulting from this procedure (Figures 15-17) showed almost 
no change compared to enhanced climatology computed from 
unfiltered data. The mean and RMS errors (Tables 5, 7) also 
show no Significant reduction. The resus wei Bowe cme Gree 
indicate that filtering the data will reduce the errors of 
simple enhancement. Filtering, however, will produce 
temperature profiles which represent a larger domain than a 
profile at a single point. Although the problem of extrapo- 
lating MLD to a new location is still unresolved, the 
concept of a "representative" MLD might preclude the require- 
ment of a MED at a single point. "“@esting Gimearc idea, 
however, will not be addressed in this thesis. 

The results of this and preceding experiments suggest 
that it is feasible to extrapolate the observed temperatures 
of most of the water column to a new location. The problem 
of determining the MLD should be dealt with separately. 
Since these experiments are relevant only to a specific time 
and domain, additional studies under different conditions 
are necessary to assess properly the usefulness of simple 


enhancement. 


48 


Table 7. Mean and RMS temperature error of filtered 
observations minus enhanced climatology computed from 
filtered observations (see text for explanation). 


Mean Error °C With RMS Error °C With 
W564) CONL1 cence S52 Confit cence 
Intervals Intervals 


EEPIn: O m 

EONS 3-56 me A—. 17 ,—. 56) 774563, 504) 
POINTS 33-68 eco, 10, —- Ca) ea Os) 
FOmms 61-96 se Z20s24,.. 41) p40. 32,049) 





eeeO/ (—2 Ay = 1 Zane Gls oye...) 3) 
ple (=. 427 ee 0) Oe (2d eee 3) 


e27(cil,.44) BD Oe) 





DEPTH: 100 m 
POINTS 3-38 ~.17(-.22,-.11) a2 Ope SO) 
POINTS 33-68 sa 3o4 427) w561( 543,270) 


POINTS 61-96 woos. 548) .51(.39,.63) 





Uae, 202) Gn 121) 


oH 2 lees 5 Oi =a 25) gO (a 32 ) 


aS 2 265— 203) D ite26 7.49) 
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IV. CONCLUSIONS 


A method of extrapolating an observed temperature profile 
from one location to another was investigated. The technique, 
referred to as simple enhancement, was examined as a possible 
application to real-time Navy ASW operations. The extrapola- 
tion required the use of a trial value which was obtained 
from two sources: an ocean thermal climatology and a real- 
time ocean forecast model (TOPS). An enhanced climatological 
temperature profile was obtained by adding an observed anomaly 
ec. , Obser ation minus climatology) to the climatology at 
some other desired location. The same method was used to 
obtain an enhanced forecast temperature profile. This proce- 
dure was shown to be a special case of optimum interpolation 
when only one observation is used, the noise-to-signal ratio 
1s zero, and the autocorrelation function between the two 
points in question (i.e., the point of the observation and the 
point where the enhanced profile is desired) is one. 

The feasibility of simple enhancement was evaluated by 
calculating mean and RMS errors of enhanced and unenhanced 
error fields. Vertical transects of the error fields were 
contoured to show the magnitude and the distribution of errors 
associated with simple enhancement. Vertical profiles were 


also ex onnenece 


2). 


The results were composited by averaging the mean and 
RMS errors for each section and depth. The average mean 
error for enhanced climatological and TOPS errors was nearly 
an order of magnitude smaller than the average for unenhanced 
climatology. The average RMS errors for the enhanced fields 
were about 1/2 of that for unenhanced climatology. 

Examination of the contoured error fields showed the 
Magnitude of error to be 1 to 2°C at the base of the mixed 
layer. In addition, temperature inversions were introduced 
in some of the enhanced temperature profiles. This is an 
indication that the MLD could not be properly extrapolated. 

The choice of the trial value did not appear critical. 
Enhanced vertical temperature profiles of climatology and 
TOPS appeared nearly identical, even though the unenhanced 
values were different. Simple enhancement was also tested 
using known spatial statistics to compute an observational 
weight from the Optimum Interpolation formalism. The errors 
in the resulting OI enhanced climatology did not show a 
reduction compared to simple enhanced fields. This is 
probably explained by the use of a single observation to 
compute the enhanced temperature profile. The OI technique 
guarantees minimum error only in an ensemble-mean sense uSing 
many observations. Thus, the approach of incorporating spatial 
statistics might not be appropriate when a single observation 


is used. 
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Finally, filtering the data prior to enhancement provides 
a method of producing a more representative temperature 
profile for a region. Enhanced climatology produced from 
filtered data shows no reduction in error over enhanced clima- 
tology produced from unfiltered data, but shows less error than 
unenhanced climatology. Thus, it is possible to produce an 
enhanced temperature profile that is also representative of a 
region, rather than a single point. The disadvantage of 
filtering is that more than a single observation is required. 

This method could provide the Navy with an improved ASW 
capability. Screening of the results, however, is recommended 
and suggests the development of an objective technique for 
removing erroneous temperature inversions. Variations of 
Simple enhancement are possible. In the present scheme, 
temperatures at each depth are extrapolated forward in space. 
A temperature profile, however, can be defined by parameters 
Other than temperature and depth pairs. A reasonable varia- 
tion would be to define a temperature profile by SST, MLD, 
and temperature gradients below the MLD. More parameters 
could be defined and added if further detail in the tempera- 
ture profile is desired. Then the observed SST, MLD and 
temperature gradients could be used to compute anomalies of 
these parameters and added to trial values of the same param- 
eters at a new location. An enhanced temperature profile 
reconstructed from these new parameters would not show temper- 


aclice —lmeierslons. 


3). 


An adequate climatology is all that is necessary for 
simple enhancement. The use or TOPS forecast Pret es (eas 
equally well, but a better use of TOPS might be to predict 
changes in MLD and SST. This would reduce eree mn tei 
time scale changes by knowing if conditions will persist or 


change in time. 
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APPENDIX A 


PRGGERING THE DATA 


Consider a small scale fluctuation at depth Z (Figure 18) 
displaced from the mean depth, Za! by afdistance DZ. ise 
therms at depth Z0 above and below Z oscillate about their 
own mean positions ZO: with departures, DZ» that decay 
toward zero at the surface and some depth below Z. The 


effects of the fluctuations can be removed by finding 2 the 


Oi” 
depth of each isotherm and replacing it with the mean depth, 


200) cf that isotherm. Assume the departures from Zn are 


proportional to departures from Z20_,- Then DZ. and DZ, are 
bevaced bys 
Dees 2 G2) Dae ey) 


Pewee paOprlace waa) 1S: 


ae) 


(D-Z,) /(D-Z_ ) La (12) 


F(Z) 


Z0/Zm Z0<2m (13) 


where D is the depth to which the effects decay to zero. 


Notice that F(Z) goes to zero as Zg approaches D, or as 


ioe aedeteomenesouriace. By detanition, 
DZo = Z0m = Zo ’ and (14) 
DZ. = Z = @. (iS) 
m m 


a7 





Figure 18. Sinusoidal representation of small scale 
fluctuations at base of mixed layer. (See text for 
explanation of variables.) 
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Then: 


ZQ,, = 2g + (29/2) (2m - 34) Zqg <2, (lke) 
Coreen ) (D-2q) (2 eS Z | Bee 2, (17) 
This formal approach is easily applied. First, select an iso- 
therm at the base of the mixed layer. Determine the mean 


depen Z4. over five adjacent temperature profiles using a 
five point running mean horizontally. Next find Z, the depth 
of the selected isotherm for each profile. Then for each 
evel 4G: compute the corrected depth 20m using equation (16) 
Cent) se teow ecoWletnGg Orotrle will nowShavem@mew corrected 
depths at each level. The temperatures remain unchanged and 


Only their vertical position has been altered. 


5) | 
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